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ABSTRACT 

We derive the effects of dark matter (DM) decays and annihilations on structure forma- 
tion. We consider moderately massive DM particles (sterile neutrinos and light DM), 
as they are expected to give the maximum contribution to heating and reionization. 
The energy injection from DM decays and annihilations produces both an enhance- 
ment in the abundance of coolants (H2 and HD) and an increase of gas temperature. 
We find that for all the considered DM models the critical halo mass for collapse, m cr i t , 
is generally higher than in the unperturbed case. However, the variation of m cr it is 
small. In the most extreme cases, i.e. considering light DM annihilations (decays) and 
halos virializing at redshift z„ r > 30 (z vir ~ 10), m crit increases by a factor ~ 4 (~ 2). 
In the case of annihilations the variations of m cr u are also sensitive to the assumed 
profile of the DM halo. Furthermore, we note that the fraction of gas which is retained 
inside the halo can be substantially reduced (to w 40 per cent of the cosmic value), 
especially in the smallest halos, as a consequence of the energy injection by DM decays 
and annihilations. 
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1 INTRODUCTION 

One of the fundamental questions concerning the formation 
of first structures is the minimum halo mass (critical mass, 
Tn c rit) for collapse at a given redshift (Silk 1977; Rees & 
Ostriker 1977; White & Rees 1978; Couchman 1985; Couch- 
man & Rees 1986; de Araujo & Opher 1988, 1991; Haiman, 
Thoul & Loeb 1996). 

Tegmark et al. (1997; T97) thoroughly addressed such 
question, pointing out how m cr it crucially depends on the 
abundance of H2, the main coolant present in the metal 
free Universe. Subsequent studies (Abel et al. 1998; Fuller & 
Couchman 2000; Galli & Palla 1998, 2002; Ripamonti 2006) 
refined the model of T97, accounting also for minor effects, 
such as the cooling induced by HD molecules. 

The production of molecules and m cr i t are sensitive to 
any physical process which can release energy in the inter- 
galactic medium (IGM). In fact, the injection of energy in 
the IGM can either delay the collapse of first halos (because 
of the increased gas temperature, or of photodissociation of 
molecules) or favour structure formation (because of the en- 
hancement in the abundance of free electrons, which act as 
catalysts for the formation of molecules). 

For this reason, it is crucial to understand the influence 
of reionization sources on structure formation. Many stud- 
ies have shown that massive metal free stars are efficient in 



dissociating H2 molecules, quenching star formation in the 
first halos (Haiman, Rees & Loeb 1997; Ciardi, Ferrara & 
Abel 2000; Ciardi et al. 2000; Haiman, Abel & Rees 2000; 
Ricotti, Gnedin & Shull 2002; Yoshida et al. 2003). Inter- 
mediate mass black holes, produced by the collapse of first 
stars, are thought to efficiently re-heat the IGM, increasing 
m cr it and reducing star formation in the smaller mass halos 
(Ricotti & Ostriker 2004; Ricotti, Ostriker & Gnedin 2005; 
Zaroubi et al. 2006). 

Also particle decays and annihilations can be sources of 
partial reionization and heating (see Mapelli, Ferrara & Pier- 
paoli 2006 and references therein) , and could influence struc- 
ture formation. For example, Shchekinov & Vasiliev (2004) 
investigated the possible effect on m cr u due to ultra-high en- 
ergy cosmic rays (UHECRs) emitted by particles decaying in 
the early Universe. Biermann & Kusenko (2006) considered 
the impact on structure formation due to sterile neutrino 
decays. Both these studies found a substantial enhancement 
on the abundance of molecular coolants (H2 and/or HD). 
However, they neglected the possible increase of gas tem- 
perature due to UHECRs or decays, respectively. 

More recently, Stasielak, Biermann & Kusenko (2006) 
evaluated the effect of sterile neutrino decays accounting also 
for the heating of the gas. However, their single-zone model 
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is likely to oversimplify the crucial behaviour of gas density 
during the halo collapse. 

In this paper, we consider the influence of dark matter 
(DM) decays and annihilations on structure formation, tak- 
ing into account variations induced both in the chemical and 
in the thermal evolution of the IGM and of the gas inside 
halos. Furthermore, we substitute the single-zone models, 
which are commonly adopted in previous papers (Haiman 
et al. 1996 is an important exception), with more sophisti- 
cated 1-D simulations. We focus on relatively low mass DM 
particles, such as sterile neutrinos and light DM (LDM), as 
their effect on the IGM is expected to be much more im- 
portant than that of heavier ( £ 100 MeV) DM particles 
(Mapelli et al. 2006). 

Sterile neutrinos are expected to decay into active neu- 
trinos and keV- photons (Dolgov 2002), while LDM can ei- 
ther decay or annihilate producing electron-positron pairs 
(Boehm et al. 2004; Hooper & Wang 2004; Picciotto & 
Pospelov 2005; Ascasibar et al. 2006). keV- photons in- 
teract with the IGM both via Compton scattering and 
photo-ionization; instead, the electron-positron pairs un- 
dergo inverse Compton scattering, collisional ionizations, 
and positron annihilations (Zdziarski & Svensson 1989; 
Chen & Kamionkowski 2004; Ripamonti, Mapelli & Ferrara 
2006, hereafter RMF06). RMF06 derived the fraction f abs {z) 
of energy emitted by sterile neutrino decays and LDM de- 
cays or annihilations which is effectively absorbed by the 
IGM through these processes. In this paper, we adopt the 
fits of fabs(z) given by RMF06. 

In Section 2 we describe the hydro-dynamical code used 
to derive and the DM models which we adopt. In Sec- 

tion 3 we discuss the effect of DM decays and annihilations 
on the chemical and thermal evolution of the IGM, giving 
an estimate of the Jeans mass. In Section 4 we describe the 
chemical and thermal evolution of the gas inside the halos, 
deriving m cr it- In the discussion (Section 5) we address vari- 
ous points, such as the variations in the baryonic mass frac- 
tion inside the halos induced by DM decays/annihilations 
and the influence of the concentration of the DM profile. 

We adopt the best-fit cosmological parameters after the 
3-yr WMAP results (Spergel et al. 2006), i.e. Q, h = 0.042, 
fi M = 0.24, Qdm = O m - fi b = 0.198, fi A = 0.76, h = 0.73, 
H = 100 h km s" 1 Mpc" 1 . 



2 METHOD 
2.1 The code 

In order to estimate the effects of the energy injection from 
DM decays and annihilations, it is necessary to follow the 
chemical and thermal evolution of primordial gas. This can 
be done with single-zone models such as the one origi- 
nally described in T97, or the adaptations by Shchekinov & 
Vasiliev (2004), and by Stasielak et al. (2006). However, this 
kind of models is bound to use some approximations which 
can be very crude. First of all, single-zone models cannot 
follow the density evolution of a virialized halo because of 
their lack of "resolution", so that it is usually assumed that 
after virialization the gas density in a halo is both uniform 
in space and constant in time. Even in the linear phase of 
the collapse, the top-hat model provides a reasonable de- 
scription of the DM component, but it becomes a rough 



approximation when the baryonic component is considered, 
as hydro-dynamical effects are likely to become important at 
scales below the Jeans or the filtering length (Peebles 1993; 
Gnedin 2000). 

An exact description would require a 3-D simulation 
(such as those of Abel, Bryan & Norman 2002, and of 
Bromm, Coppi & Larson 2002); but it is possible to cap- 
ture the basic features of the collapse phenomenon by an 
intermediate, time-efficient approach, e.g. by means of 1-D 
simulation, as they are both more accurate than single-zone 
models, and much faster than 3-D simulations (Haiman et 
al. 1996). 

In this paper, we use the 1-D Lagrangian, spherically 
symmetric code described by Ripamonti et al. (2002), as up- 
dated in Ripamonti (2006). Such a code includes the treat- 
ment of: 

• the gravitational and hydro-dynamical evolution of the 
gas (by means of an artificial viscosity scheme) ; 

• the chemical evolution of 12 species (H, H + , H~, H2, 
H+, D, D+, HD, He, He+, He ++ , and e"; see table 1 of 
Ripamonti 2006 for a list of the considered reactions and of 
the adopted reaction rates); 

• the cooling (or heating) due to a number of components, 
such as the compressional (adiabatic) heating, the emission 
and absorption of line radiation from H, H2 and HD (ac- 
counting for the effects of the cosmic microwave background, 
CMB), the heating (or cooling) from the Compton scatter- 
ing of CMB photons off free electrons, and the heating (or 
cooling) due to chemical reactions (e.g. the formation or dis- 
sociation of H2 molecules) ; 

• the gravitational effects of DM, according to a simple 
model which is based on the top-hat formalism up to the 
turn-around redshift, smoothly evolving into a concentrated 
profile after virialization (see Section 2.2 for more details). 



2.2 DM profiles 

The code does not include a self-consistent treatment of DM. 
Instead, the function pDivi(r, z), describing the DM density 
profile and its redshift evolution, must be chosen a priori. 
Since our results might depend on this choice, we decided 
to study two quite different cases, which we call 'isothermal' 
and 'NFW (from the profile described in Navarro, Frenk & 
White 1997) depending on the shape of the density profile 
after virialization. 

In both cases the DM distribution was assumed to 
be spherically symmetric, and concentric with the simu- 
lated region, whose central part represents the halo which 
is being investigated. At any redshift a DM mass Mdm = 
-Mhaio^DM/fiM is assumed to be within the truncation ra- 
dius 



Rtr{z) 



( 3 M DM A 
t Rvi 



1/3 



t(Zvir)-'(zta) 



if Z > Z ta 

if Zta > Z > Z vir ^ 
if Z < Z v i r 



where t(z) is the time corresponding to redshift z, and z v ; r , 
Zta, -Rvir and Pth(z) are the halo virialization and turn- 
around redshifts, its virial radius, and the DM density inside 
the halo at z > z ta (as derived from the evolution of a simple 
top-hat fluctuation; see e.g. Padmanabhan 1993, or T97), 
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respectively. The exact definition of these quantities can be 
found in Ripamonti 2006. 

The isothermal and NFW assumptions differ only after 
the halo virialization (z < 2 v ; r ); both of them refer to a 
static DM profile. In the isothermal case 

{pcorc if r < R core 5 

Pc OIC (r / R cor c)~ 2 if -Rcorc < r < Rti\ (2) 
p fiD M (l + z) 3 if r>R tI ; 

instead, in the NFW case the DM density profile is chosen 
to be 



PNFW 

Po"dm(1 + z) 



if r < Rtr 
if r > R tI 



(3) 



where po — 1.88 x 10~ 29 h 2 gem -3 is the critical density of 
the Universe at present. In both cases 7? CO re = £,Rvir, where 
£ = 0.1 is a parameter (cfr. Hcrnquist 1993, and Burkert 
1995 for the choice of its value; in the NFW case £ the inverse 
of the more commonly uses concentration parameter), and 
the densities p cor c and pnfw can be found by requiring the 
DM mass within R tr to be equal to Mom- 

At z > z v i r both the isothermal and NFW case assume 
the DM density profile described by equation 2. However, 
at such redshift the profile is not static, because the core 
radius is evolved with redshift 



-^corc {z) 



Rtr(z) 



Rv 



t(z vir )-t( Zta ) 



if Z > Zta 
if Zta > Z > 



(4) 



Such a choice combines the behaviour of a top-hat fluctua- 
tion (the density inside R tr is assumed to be uniform until 
the turn-around) with a transition to the final density pro- 
files. 

We only considered the case £ = 0.1 (i.e. concentration 
10 for the NFW profile), instead of varying £, because the 
differences between the isothermal and NFW cases are quite 
relevant even with the same value of £. In fact, the NFW case 
is representative of concentrated halos, whereas the isother- 
mal case is representative of relatively shallow potentials. 

However, it is important to note that the flat central 
profile of the isothermal case helps to ensure that the be- 
haviour we observe near the centre is due to the self-gravity 
and hydrodynamics of the simulated gas, rather than to the 
assumed DM profile. 



2.3 Treatment of the DM energy injection 

The above code was modified in order to include the effects 
of the energy injection from DM decays/annihilations on 
both the chemical and thermal evolution of the gas. 

The gas can be heated, excited and ionized by the en- 
ergy input due to DM decays/annihilations. It is important 
to note that the fraction of the absorbed energy going into 
each one of these components is quite unrelated to how the 
energy was deposited in the IGM at the first step. For exam- 
ple, if a keV a photon ionizes an atom, the resulting electron 
will generate a cascade of collisions, and the energy of the 
photon will go not only into ionizations, but also into exci- 
tations and heating. 

Thus, given the energy injection per baryon from DM 



Table 1. List of chemical reactions stimulated by the DM energy 
injection. 



Species 


Reaction 


Threshold 


H 


H -> H+ + e" 


13.6 eV 


He 


He -> He+ + e" 


24.6 cV 


He+ 


He+ -> He++ + e" 


54.4 eV 


D 


D -> D+ + e" 


13.6 eV 


H 2 


H 2 ->• H + H 


4.48 cV 


HD 


HD -> H + D 


4.51 eV 


H 2 + 


H+ -» H + H+ 


2.65 eV 



decays and annihilations, e(z) (described in Section 2.4), we 
split it into an heating and ionization component 1 



Cheat (z) — 
eion(z) = 



C[l~(l-x(z)f] b e(z) 
1 - x(z) 

n < Z )> 



(5) 
(6) 



where x(z) is the ionization fraction. In the first equation 
we are using the fit to the results of Shull & Van Steenberg 
(1985) which is provided in their paper (with C = 0.9971, 
a = 0.2663, and b — 1.3163), while in the second equation 
we are using the fit to the same results given by Chen & 
Kamionkowski (2004). 

The heating component is simply added to the equa- 
tions describing the thermal state of the gas. Instead, the 
ionization component is further split between H, He, He + , 
D, H 2 , HD and Hj , according to their number abundance: 



*l2je{H,He,He+ ,D,H 2 ,HD,H+) 



.(*) 



(7) 



where the indices i and j indicate chemical species, and JVj = 
Ni if the species i is atomic, or as ATj = 2iV» if the species i is 
molecular; Ni is the number density (per unit volume) of the 
chemical species i. In principle, the terms in the sum above 
should be weighted by the cross-section for each species. 
However both the cross-section and the energy spectrum are 
too complex to be implemented in our simple calculations. 
In particular, the energy spectrum is expected to be the 
result of a cascade (Shull & Van Steenberg 1985). The error 
in neglecting these factors is quite small, as m cr it is more 
sensitive to the temperature increase than to the chemistry 
(see Section 4). 

The quantity e; approximates the energy which is ab- 
sorbed by chemical reactions dissociating the species i, and 
is translated into a reaction rate (number of reactions per 
particle per unit time) through a division by the energy 
threshold Eth,i of the considered reaction. The list of the 
reactions and of the energy thresholds is given in Table 1. 

We neglect the absorption of ionization energy by H + , 
He ++ , D + and e~ , because these species cannot be ionized 
further. We also neglect the ionization energy absorption by 
H~ because the energy threshold for the transformation of 

1 The sum of these two components is less than e(z), as a sig- 
nificant fraction of the injected energy goes in atomic/molecular 
excitations and does not affect the chemical or thermal state of 
the gas. 
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H~ into H is negative, and cannot be treated with our simple 
formalism. However the number abundance of BP is always 
very small, and the number of dissociations induced by DM 
decays/annihilations is likely to be negligible. 



2.4 DM models 

We apply this formalism to two different DM candidates, 
i.e. sterile neutrinos and LDM. Sterile neutrinos are one of 
the most popular warm DM (WDM) candidates (Colombi, 
Dodelson & Widrow 1996; Sommer-Larsen & Dolgov 2001; 
Bode, Ostriker & Turok 2001). They can decay via different 
channels (Dolgov 2002; Hansen & Haiman 2004). In this 
paper we are interested in the radiative decay, i.e. the decay 
of a sterile neutrino into a photon and an active neutrino, 
because of its direct impact on the IGM (Mapelli & Ferrara 
2005; Mapelli et al. 2006). The photon produced in the decay 
interacts with the IGM both via Compton scattering and 
photo- ionization (RMF06). 

LDM particles have recently become of interest, be- 
cause they provide a viable explanation for the detected 
511-keV excess from the Galactic centre (Boehm et al. 2004; 
Knodlseder et al. 2005). If they are source of the 511-keV 
excess, then their maximum allowed mass ttildm should be 
20 MeV, not to overproduce detectable gamma rays via in- 
ternal Bremsstrahlung (Beacom, Bell & Bertone 2004). If 
we consider also the production of gamma rays for inflight 
annihilation of the positrons, this upper limit might become 
~ 3 MeV (Beacom & Yiiksel 2006). 

LDM can either decay or annihilate, producing photons, 
neutrinos and pairs. In this paper we consider both LDM 
decays and annihilations, but we restrict their treatment to 
the case where the product particles are e + — e~ pairs. In 
fact, in the case of pair production the impact of LDM on the 
IGM is maximum (RMF06). The e + — e~ pairs are expected 
to interact with the IGM via inverse Compton scattering, 
collisional ionization and positron annihilation (RMF06). 



2.4-1 The energy input from the "background" 

We first consider the energy injected in the general IGM 
after cosmic DM decays/annihilations. 

Both in the case of sterile neutrinos and of LDM, the 
rate of energy transfer (per baryon) to the IGM because of 
this "background" contribution can be written as: 



£bkg(z) = fabs(z) HDM^) TTlOM C , 



(8) 



where itidh is the mass of a DM particle and c is the speed 
of light. The energy absorbed fraction, f a bs{z), has been de- 
rived in RMF06; wdm(z) is the decrease rate of the number 
of DM particles per baryon. 

In the case of DM decays, ?idm(z) is given by 



nuM(z) 



riDM,0 
TDM 



(9) 



where n D M,o and t D m are the current number of DM par- 
ticles per baryon and the lifetime of DM particles, respec- 
tively, t dm is assumed to be much longer than the present 
value of the Hubble time, as is the case for all the models 
we are considering. 

For the annihilations: 



huM(z) ~ -n DMfi N b (0){av)(l + z) 3 , (10) 

where Nt(0) = 2.5 x 10 -7 g cm -3 is the current baryon num- 
ber density (Spergel et al. 2006), and (o~v) is the thermally 
averaged DM annihilation cross-section. 

Both for sterile neutrinos and LDM, the values of noM.o, 
tdm and (av) adopted in this paper are the same reported 
in RMF06; for convenience, they are listed in Table 2. 



2.4-2 The "local" energy input 

In addition to the background energy injection discussed 
above, the baryons inside a halo absorb extra energy from 
the additional decays/annihilations of overdense halo DM. 
In the case of decays the total excess energy "produced" 
inside the halo is 

47riV b (0)(l + 2) 3 mDMC 2 



Ei oc (z) = 



TDM 



R tr (z) 



drr [nDM(r, z) — tidm,o], 



(11) 







and in the case of annihilations 



Ei oc (z) 



4tt [N b (0)(l + zf] 2 m DM c (av) 



I 



dr-r 2 [n DM (r,z) 



«dm,o(z)L (12) 



where n D M(r, z) = pDu(r, z) /[m OM N b (0)(l + z) 3 ] is the 
number of DM particles per baryon at redshift z and at 
a distance r from the centre of the halo. 

Then, we compute the baryon column density Ej, from 
the halo centre to the truncation radius, and find the frac- 
tion of the energy Ei oc which is absorbed by such a column 
density, /; oc . If the DM produces photons of energy E-y, 

f loc = E 6 \a He+H (E^ + * T -%g f-^j)l (13) 

m e c 2 \m e c 2 / . 

where OHe+H{E) is the photo- ionization cross section (see 
Zdziarski & Svensson 1989; RMF06), o T is the Thomson 
cross section, and the function g is defined in equation 4.9 
of Zdziarski & Svensson (1989). Instead, if the DM produces 
electron-positron pairs with Lorentz factor 7 (and energy 
E=7 m e c 2 ), 



flo 



rg (f>e,ion(z, E) + 



Rtr 



-<j> e , com (z,E)(U) 



cJV 6 (0)(l 

where 4> e ,ion(z, E) and 4>e,c°m(z, E) give the fraction of the 
energy of an electron/positron which is absorbed per unit 
time by the IGM, because of collisional ionizations and 
of Compton scattering of CMB photons, respectively. Such 
functions are given by equations 14-15, and 16-18 of RMF06. 

We then assume that all the baryons inside the halo ab- 
sorb the same amount of energy from this local contribution. 
So, the "local" energy deposition in each baryon within the 
truncation radius Rt r (z) is 



£ioc(z, r) = Elocfloc 



M gas [Ru(z)] 



(15) 



where M ga s(Rtr) is the mass of the gas inside the trunca- 
tion radius, and ran is the mass of an H atom. Instead, for 
baryons at a distance larger than R tr (z) from the centre of 
the halo we assume ei oc (z,r) = 0. 
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Table 2. List of the main characteristics of the considered DM models. In particular, from the leftmost to the rightmost column: DM 
model, 71dm,0: T DM (if decaying particle), (crv) (if annihilating), comoving free-streaming lengths (Aps jn and Aps,i) and the associated 
mass scale. Aps is taken to be the maximum between Apg jn and Apg^. 



DM model 


«DM,o/10 3 


tdm/(10 27 s) 


(at>)/(10- 29 cm 3 s- 1 ) 


AFS.n(pc) 


AFS,i(P c ) 


M FS (M ) 


v 4 keV (decaying) 


1176 


2.96 




3 x 10 5 


60 


5 x 10 s 


^15 keV (decaying) 


314 


1.98 




8 x 10 4 


35 


1 x 10 7 


v 25 keV (decaying) 


188 


0.097 




4.8 x 10 4 


25 


2 x 10 6 


LDM 3 MeV (decaying) 


1.49 


1.2 




2.2 


< 0.1 


2 x 10~ 7 


LDM 10 MeV (decaying) 


0.446 


4.0 




0.5 


< 0.1 


3 x 10~ 8 


LDM 1 MeV (annihilating) 


4.46 




4 


9.5 


300 


0.5 


LDM 3 MeV (annihilating) 


1.49 




12 


2.2 


100 


0.02 


LDM 10 MeV (annihilating) 


0.446 




24 


0.5 


30 


5xl0~ 4 



The total energy input per baryon from DM de- 
cays/annihilations is then 

e(r, z) = e bkg (z) + e ioc (r, z). (16) 



2.5 The free-streaming/damping lengths 

We need to take into account that small DM fluctuations 
might be washed out by damping mechanisms, and in 
particular by free-streaming (see e.g. Padmanabhan 1993; 
Sommer-Larsen & Dolgov 2001; Boehm et al. 2005). 

The free-streaming scale depends on the considered par- 
ticle, and on the strength of its interactions. If such inter- 
actions are negligible, the comoving free-streaming length 
is 

f 0.3 ~' ( <PC/ 3 ( if T>> ) Mpc for neutrinos 

-^FS,n ^ < , . -4/3 

[ 0.20(fi DM ft 2 ) 1/3 (=j»gv-J Mpc for LDM, 

where ttidm is the mass of the considered DM particle, p and 
T are the modulus of the proper momentum and the tem- 
perature of the particle, respectively. The above expressions 
are derived from Abazajian, Fuller & Patel (2001), and from 
Boehm et al. (2005), respectively. 

However, if the DM particle interacts at a non negligi- 
ble rate, the comoving free-streaming length is (Boehm et 
al. 2005; Boehm & Schaeffer 2005) 

with 

-L dec, DM = TdM (Zdec)(l+Zdec) , (19) 

where r D ivi(.Zdcc) is the DM interaction rate at its decoupling 
redshift Zdcc- 

In the case of decaying particles, we assume that this 
rate is simply the inverse of the present lifetimes of DM 
particles; so, rDM(^dec) < 10~ 24 s _1 in all the cases we are 
considering (cfr. Table 2). Since z dcc S> 10 3 , equations (17) 
and (18) imply that Aps,i <C Aps, n ; so, for decaying particles 
we will use a free-streaming scale Afs = AFs, n . 

Instead, when annihilating particles are considered, the 
interaction rate is r D M(zdoc) = \ nDM,o-Nb(0)(<7u) dcc (l + 
Zdec) 3 , where (av)dcc ~ 10~ 26 cm 3 s^ 1 (see, for example, 



Ascasibar et al. 2006) is the thermally averaged annihila- 
tion cross section at the epoch of decoupling 2 . Then, in the 
case of annihilating particles, AFs,i 3> Aps, n , and we must 
assume Afs = Afs,i- In Table 2 we give the detailed list of 
the free-streaming lengths for each DM model. The same 
table also lists the free-streaming mass scale 

Mfs = t(^t) 3 poflMh2 - (20) 

Damping erases fluctuations on scales smaller than Aps- 
Thus, objects of mass Mfs are unlikely to form, unless 
the small scale power spectrum is regenerated by non-linear 
effects at low redshift (see Boehm et al. 2005; however, this 
regeneration appears to take plate at z ~ 2, which is much 
later than the epoch we are considering). To introduce cor- 
rections for this effect in our models is beyond the purposes 
of this paper. 

Damping might also affect the density profile of halos, 
erasing cusps on scales £ Afs • The values of Afs and Mfs 
listed in Table 2 show that this effect of damping can be 
important only for sterile neutrinos. In the case of isother- 
mal density profiles, we account for it adopting the follow- 
ing correction. If the core radius R CO ic(z) (equation 4) is 
smaller than Afs/2, we increase its value to i? C oro,now(z) = 
min(Rti(z), Afs/2. Instead, in the case of simulations with 
NFW density profile we do not introduce this correction, be- 
cause these simulations are intended to explore the effects 
of high concentration (see Section 4.2). 

2.6 The simulations 

Our code was used to run a large number of simulations, in 
order to explore a wide range of the z vir — M ha i parameter 
space: we considered halo masses in the range 10 4 -10 7 M Q , 
and virialization redshifts between 10 and 100. 

We actually simulate a mass which is 1000 times higher 
than that of the collapsing halo, in order to include in the 

2 For annihilating particles, a third scale length might be con- 
sidered, i.e. the mixed damping length \ m< i (Boehm et al. 2003). 
The upper limit of A m( j is 6.7, 2.9 and 1.2 kpc for 1, 3 and 10 
MeV LDM, respectively. These values are a factor ~ 20 larger 
than the corresponding Aps,i; but our conclusions about m cr n 
(sec Section 4) are not significantly affected, because Mps is still 
smaller than the mass range we are considering. 
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Redshift 

Figure 1. Effects of decaying/annihilating DM on the IGM evo- 
lution. Left axis: fractional abundances of free electrons (e — ), H2 
and HD as a function of redshift. Right axis: matter temperature 
as a function of redshift. Top panel: Effect of decaying sterile neu- 
trinos of mass 25 keV (dashed line). Central panel: decaying LDM 
of mass 10 MeV (dashed line). Bottom panel: annihilating LDM 
of mass 1 MeV (dashed line). In all the panels the dotted line is 
the CMB temperature and the solid line represents the thermal 
and chemical evolution without DM decays/annihilations. 



simulation a mass which is larger than the cosmological fil- 
tering mass (Gnedin 2000); otherwise, our treatment of hy- 
drodynamics might be incorrect. The simulated object is 
divided in 150 shells, whose spacing was chosen so that (i) 
the mass of the shells smoothly increases when moving out- 
wards, (ii) the central shell always encloses a gas mass of 
~ 0.3Mq, and (iii) the central 100 shells initially enclose a 
mass ~ Mhaio (including the DM). 

The simulations are started at z — 1000, when we as- 
sume that the gas density is uniformly equal to the cosmo- 
logical value (whereas the DM density profile is not perfectly 
uniform; see Section 2.2), the gas temperature is equal to 
the CMB temperature (~ 2.728 K), and adopt the chemical 
abundances listed in table 2 of Ripamonti 2006 3 . 

They are stopped either when the gas density reaches 
the threshold pjmu = 10 5 cm -3 , or after a time 2t(z v i r ) has 
elapsed. 

Each set of simulations was repeated for each different 
DM decay /annihilation model, and also for the "standard" 
case without any energy injection from DM, which is used 
as a reference against which we compare our results. 



3 The most important abundances listed there are the H ion- 
ization fraction (at z = 1000) N H+ /N H = 0.0672, the helium 
abundance N[f B /N[f = 0.0833, and the deuterium abundance 
N D /N H = 2.5 x 10- 5 . 



3 IGM EVOLUTION 

In order to look at the influence of DM decays and annihi- 
lations upon the IGM, we have used a simplified version of 
our code (where the density was assumed to evolve as the 
cosmological value). 

3.1 Chemistry and temperature 

Our code follows the chemical evolution of 12 chemical 
species (see previous Section). Two of them, molecular hy- 
drogen (H2) and HD are particularly important for our pur- 
poses, because they are the main coolant of the metal free 
gas. In Fig. 1 we show the fractional abundances of both H2 
and HD together with the ionized fraction and the matter 
temperature as a function of redshift for the considered DM 
models. 

In all the DM models, both the matter temperature and 
the abundance of H2, HD and free electrons are enhanced 
by DM decays/annihilations. This effect is smaller for sterile 
neutrinos than for LDM particles. The main difference be- 
tween decays and annihilations is represented by the redshift 
range in which the influence of DM is important. For all the 
considered quantities (i.e. temperature and abundance of 
e~, H2 and HD) the energy injection from DM decays starts 
to be significant at redshift lower than ~ 100. Instead, the 
influence of annihilations is important already at redshift 
~ 900. The annihilations represent also the case where the 
abundance of the two coolants is most enhanced (a factor 
~ 17 for the H2 and ~ 90 for the HD). Furthermore, the an- 
nihilations keep the matter temperature close to the CMB 
temperature everywhere, up to z ~ 50 — 100. This fact can 
have important consequences for experiments searching for 
high redshift HI 21-cm line signals (Shchekinov & Vasiliev 
2006). 

3.2 Jeans mass 

In order to establish the influence of DM decays or anni- 
hilations on the structure formation, the key point is the 
following. DM decays and annihilations increase both the 
matter temperature and the abundance of coolants. The for- 
mer effect tends to delay the formation of structures, while 
the latter favours an early collapse of the halos. Which of 
these two opposite effects is dominant? When looking at the 
average properties of the IGM, the most popular diagnostic 
is the cosmological Jeans mass, m.j (Peebles 1993): 

frri \ 7T / TYk B T \ i/2 _ 1/2 

mj(T ' P ' M) = 6 {g^J 9 

~ 5OAf T 3/ V 3/2 (— ) , (21) 
\rriH J 

where ks is the Boltzmann constant, G the gravitational 
constant and fi the mean molecular weight. 

In Fig. 2 we show the evolution of mj. For all the consid- 
ered DM models mj is considerably increased by the effect 
of decays and annihilations. This means that, when the IGM 
is considered, the increase of the matter temperature dom- 
inates over the enhancement of the coolant abundance. For 
example, at 2 = 10 nij = 3.0 x 1O 7 M0 in the case of 10-MeV 
LDM decays, a factor ~ 420 higher than in the unperturbed 
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Figure 2. Jeans mass as a function of redshift. Top panel: Effect 
of decaying sterile neutrinos of mass 25 (dashed line), 15 (dot- 
ted) and 4 keV (dot-dashed). Central panel: Decaying LDM of 
mass 10 (dashed line) and 3 MeV (dotted). Bottom panel: An- 
nihilating LDM of mass 10 (dashed line), 3 (dotted) and 1 MeV 
(dot-dashed). In all the panels the solid line represents the ther- 
mal and chemical evolution without DM decays/annihilations. 
The horizontal dot-dashed line shows the behaviour of the Jeans 
mass calculated by assuming that the gas temperature is always 
equal to the CMB temperature. 



case (for which mj = 7.2 x 10 4 Mq). The increase of mj is 
less pronounced, but nevertheless significant in the case of 
sterile neutrino decays (mj = 6.6 x 10 s Mq for 25-keV ster- 
ile neutrinos) and LDM annihilations (mj = 3.9 x 10 5 Mq 
for 1-MeV LDM). 

From this fact one could naively infer that DM decays 
and annihilations strongly delay the formation of the first 
structures. However, m.j refers only to the global proper- 
ties of the IGM, and does not account for the non-linear 
evolution of collapsing halos. 



4 EVOLUTION INSIDE HALOS 
4.1 Chemistry and temperature 

To assess the star forming ability of the first halos, it is 
necessary to follow their hydro-dynamical and chemical evo- 
lution. This has been done by using the code described in 
Section 2. Fig. 3 (Fig. 4) shows, as an example, the evolution 
of the most relevant properties of the gas at the centre of a 
6 x 10 5 Mq (2 x 10 6 Mq) halo, virializing at z vir = 12. In 
these figures the effects of different DM decay /annihilation 
scenarios are compared with the unperturbed case. For the 
smaller halo (Fig. 3) LDM annihilations and sterile neutrino 
decays delay the collapse by ~ 1 redshift units, and in the 
case of LDM decays this effect is even more significant: in 
presence of 10 MeV-LDM decays the halo has not collapsed 
yet after one Hubble time from virialization. Instead, for the 
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Figure 3. Evolution of the central region of a 6 x 10 Mq isother- 
mal halo virializing at z v i r = 12. From left to right and from 
top to bottom: density, temperature, electron abundance and H2 
abundance as function of redshift. The solid line represents the 
unperturbed case (i.e. without DM decays and annihilations). The 
dashed, dot-dashed and dotted lines account for the contribution 
of 25-keV sterile neutrino decays, 1-MeV LDM annihilations and 
10-MeV LDM decays, respectively. 
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Figure 4. Evolution of the central region of a 2 X 10 Mq isother- 
mal halo virializing at z v i r = 12. From left to right and from 
top to bottom: density, temperature, electron abundance and H2 
abundance as function of redshift. The dashed, dot-dashed and 
dotted lines are the same as in Fig. 3. 
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larger halo (Fig. 4) the difference between the case with and 
without DM decays/annihilations is negligible. 

It is worth to discuss in detail the gas temperature (T) 
and number density (TV = p/ran) evolution shown in Fig. 3. 
At the epoch of turn-around (z ~ 19 — 20) all the models 
have the same TV; but in the unperturbed one T is 1.5-5 
times lower than in the others. This is important, because 
the fast increase in T and TV due to the virialization pro- 
cess is essentially adiabatic, as it happens on a time-scale 
much faster than that for cooling. In the adiabatic approxi- 
mation, r ad = -PdV oc PIN' 1 ~ (TV + ATV) -1 ] (where P 
is the gas pressure), is proportional to the initial T for any 
given increase in TV. T and P increase faster in models "pre- 
heated" by DM decays/annihilations, so that the pressure 
gradient slowing and halting the collapse develops earlier. 
In fact, the phase of unimpeded collapse stops at TV ~ 10 
cm" 3 for the unperturbed case, but only at TV ~ 0.3 — 3 
cm -3 for the other models. At these densities the cooling 
per unit mass is proportional both to the H2 fraction and 
to TV, so that the higher central density of the model with- 
out DM decays/annihilations largely compensates its lower 
H2 abundance, i.e. the unperturbed case is the fastest cool- 
ing one. However, it is interesting to note that the final T is 
much lower for the cases where DM energy input is included, 
because of the enhancement of the HD fraction, which pro- 
vides an efficient cooling mechanism also at T ,$ 200 K (Ri- 
pamonti 2006). 

4.2 Critical mass 

Each of our simulated halos was classified as collapsing 
(or, equivalently, efficiently cooling) or non-collapsing (in- 
efficiently cooling), depending on whether it reaches a max- 
imum density larger than 10 5 mu cm -3 in less than a Hub- 
ble time. We define the critical mass, m C ru, as the minimum 
mass of a collapsing halo at a given virialization redshift z V i r - 

The values of m cr it as a function of the virialization 
redshift are shown in Fig. 5 for isothermal halos, and in 
Fig. 6 for NFW halos. The main trend inferred from the 
behaviour of mj (i.e. the delay of structure formation) is 
confirmed by merit- However, the influence of DM decays 
and annihilations on m cr it is much smaller than could be 
expected from mj. 

In the case of isothermal halos, the effects of the DM 
decays on m cr it (top and central panel) are completely negli- 
gible at high redshift and become significant only at redshift 
less than 20. Even for z V i r = 10 and 10-MeV decaying LDM 
(which has the strongest effect), m cr it increases only of by 
a factor 2. 

Sterile neutrinos (especially with mass of 4 keV) seem to 
significantly increase merit- However, this is not connected 
with sterile neutrinos decays; but it is mainly due to dif- 
ferences in the assumed DM density profile, as the free- 
streaming length for these WDM halos is larger than the 
"standard" value of the core radius (0. l_R v i r ) and the "damp- 
ing" correction described in Section 2.5 is important. In fact, 
m cr it is almost the same down to z ~ 20, if we consider 
(thick dot-dashed line in Fig. 5) or neglect (thin solid line) 
sterile neutrino decays, and the difference between these two 
cases is small even at lower redshift. 

In the case of LDM annihilations in isothermal halos, 
m cr it is higher than in the unperturbed case for every con- 
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Figure 5. The critical mass as a function of the virialization 
redshift z v i r for isothermal halos. Top panel: Effect of decaying 
sterile neutrinos. Central: Decaying LDM. Bottom: Annihilating 
LDM. The lines used in the three panels are the same as in Fig. 2. 
The thin solid line used in the top panel represents the case of 
non-decaying WDM (with mass 4 keV). 
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Figure 6. The critical mass as a function of the virialization 
redshift z„; r for NFW halos. The lines used in the three panels 
are the same as in Fig. 5. 

sidered virialization redshift (z v ir — 10 — 100), confirming 
that the annihilations play a role even at very high z. How- 
ever, the difference with respect to the unperturbed m cr i t is 
always less than a factor ~ 2. 

If we assume a NFW profile for DM halos (Fig. 6), 
the general effect is a substantial decrease of m cr i t (hide- 
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Figure 7. Halo baryonic mass fraction (see definition in equa- 
tion 22) as a function of the halo mass for a fixed virialization 
redshift (z v i r = 10) and for an isothermal DM profile. Top panel: 
Effect of decaying sterile neutrinos. Central: Decaying LDM. Bot- 
tom: Annihilating LDM. The lines used in the three panels are 
the same as in Fig. 2. Thick (thin) lines indicate that the halo 
mass is larger (smaller) than m cr jt. 

pendently of DM decays/annihilations), due to the higher 
central densities, and the consequent stronger gravitational 
pull. LDM decays tend to increase m cr i t (at z < 20) with 
respect to the unperturbed case; the increase factor is only 
marginally larger than in the case of isothermal halos. Sterile 
neutrino decays appear to slightly reduce m cr it with respect 
to the unperturbed case; but this result is unphysical, as it 
does not account for damping. It is worth to note that, even 
if we ignore damping effects and adopt the models with the 
highest possible concentration, the influence of sterile neu- 
trino decays on structure formation remains very small. 

Instead, in the case of LDM annihilations, m cr it gener- 
ally increases more in NFW halos than in isothermal halos. 
This is due to the importance of the local contribution for 
annihilations (see Section 5.2). 



5 DISCUSSION 
5.1 Gas fraction 

Why considering either mj or m crit leads to such widely 
different conclusions? As mentioned earlier, mj is a rough 
tracer of the minimum halo mass for collapse; but it is not 
sensitive to the local properties of the halo. For example, 
to calculate mj one has to assume that the density is uni- 
form. On the contrary, m cr u depends on the evolution of 
the central region of the halo, where the density increases 
much more rapidly than in the outskirts, driving the collapse 
(Figs. 3 and 4). When the density becomes sufficiently high, 
the molecular cooling largely overcomes the heating due to 
DM decays/annihilations (Figs. 3 and 4). This is the reason 
why m cr it does not increase significantly in presence of DM 
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Figure 8. Evolution of m cr it as a function of the virialization red- 
shift in the case without DM (solid lines), with 3-MeV LDM de- 
cays (dotted) and with 1-MeV LDM annihilations (dot-dashed). 
Thin (thick) lines are with (without) the local contribution of the 
DM inside the halo. Note that the thin dotted line is completely 
superimposed to the thick dotted line, because the local contribu- 
tion of decaying DM is always negligible. Top panel: isothermal 
halos. Bottom panel: NFW halos. 



decays/annihilations. Instead, mj, which does not account 
for the density increase and for the consequent cooling en- 
hancement, dramatically grows. Then, we can conclude that 
mj does not provide a realistic estimate of the effects of DM 
decays/annihilations. 

Is there any physical implication of the growth of mj? 
mj is directly related to the global hydro-dynamical be- 
haviour of the gas inside a halo: in halos with mass below 
raj the gas pressure prevents the development of a large 
gas overdensity, while in more massive halos the gas accu- 
mulation should proceed almost unimpeded. If so, the mass 
fraction in baryons within the virial radius of a halo should 
be of the order of the cosmological average Q),/Q,m when 
Mhaio ~ rrij, but it should be lower when M^aio ~ rnj. To 
check this hypothesis, we derive the baryonic mass fraction 
{f gas) normalized to Qi,/Qm, i-e. 

, _ M gas (R vir ) Sl M 

Jgas — , 
■Mhalo "6 

where M gas (R vir ) is the mass of gas within the virial 
radius at the final stage of each simulation; instead Mhato is 
the total mass of the halo at the beginning of the simulation. 

In Fig. 7 we show f gas as a function of mass inside halos 
virializing at redshift z V i r — 10. The comparison between 
the various scenarios clearly shows that the energy injection 
from DM decays/annihilations can substantially reduce the 
gas fraction inside all halos, especially the smallest ones. 
This is clearly related to the increase in mj, as the largest 
variation in f ga s occurs in the LDM decay scenario, where 
the increase of mj is maximum (Fig. 2). 
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5.2 Concentration and local contribution 

In this paper we basically considered two different DM pro- 
files: isothermal and NFW halos. These two models are com- 
parable in terms of virial radius and main properties; but the 
NFW one is much more concentrated. From the comparison 
between Fig. 5 and Fig. 6 we have already seen that this 
difference has important effects on the critical mass: m cr i t 
is generally much lower for a NFW than for an isothermal 
halo, independently of DM decays and annihilations. For 
decaying particles, higher concentrations lead to marginally 
stronger effects in delaying the formation of structures, but 
the effect is minimal. For annihilations this effect of high 
concentrations is much larger (see Figs. 5, 6 and 8). 

Another significant characteristic of our model is 
the estimate of the local contribution due to DM de- 
cays/annihilations occurring inside the halo (see Section 
2.4.2). Pointedly, the local contribution strongly depends on 
the DM profile. But what is the importance of the local con- 
tribution? In Fig. 8 we indicate the effect on m cr i t of LDM 
decays (dotted line) and annihilations (dot-dashed), by in- 
cluding (thin line) or not (thick) the local contribution. In 
the case of DM decays the thin and thick dotted lines appear 
superimposed, both in the isothermal (top panel) and NFW 
(bottom) profile, indicating that the local contribution is 
always negligible for decays. 

Instead, if we consider the annihilations, the case with 
(thin line) and without (thick) the local contribution arc 
very different, especially for the NFW profile. If we do not 
include the local contribution, m cr i t is very close to the un- 
perturbed value. If we switch on the local contribution in 
the isothermal profile, m cr i t is substantially higher (a factor 
~ 2) than in the unperturbed case, at least for high virial- 
ization redshifts (z vir > 30). Finally, if we account for the 
local contribution in a NFW halo, m cr i t is always higher 
(a factor ~ 2 — 4) than in the unperturbed case (even at 
low virialization redshifts). This is consistent with the fact 
that the annihilation rate strongly depends on the local DM 
density. 



6 CONCLUSIONS 

In this paper we derived the effects of DM decays and anni- 
hilations on structure formation. We considered only mod- 
erately massive DM particles (sterile neutrinos and LDM), 
as they are expected to give the maximum contribution to 
heating and reionization (Mapelli et al. 2006). To describe 
the interaction between the IGM and the decay /annihilation 
products we followed the recipes recently derived by RMF06. 

We accounted not only for the diffuse cosmological con- 
tribution to heating and ionization, but also for the local 
contribution due to DM decays and annihilations occurring 
in the halo itself. The local contribution results to be dom- 
inant in the case of DM annihilations especially for cuspy 
DM profiles. 

The energy injection from DM decays/annihilations 
produces both an enhancement in the abundance of coolants 
(H2 and HD) and an increase of gas temperature. We found 
that for all the considered DM models (sterile neutrino de- 
cays, LDM decays and annihilations) the critical halo mass 
for collapse, often higher than in the unperturbed 



case. This means that DM decays and annihilations tend 
to delay the formation of structures. However, the variation 
of m cr it is minimal. In the most extreme cases, i.e. consid- 
ering LDM annihilations (decays) and halos virializing at 
redshift z v i r > 30 (z v i r ~ 10), m cr i t increases of a factor 
~ 4 (~ 2). In the case of decays, the variations of m cri t are 
almost independent from the assumed concentration of the 
DM halo, although higher concentrations (corresponding to 
smaller values of m cr it) seem to be associated with slightly 
stronger effects of the DM energy injection. The dependence 
on concentration is more evident in the case of annihilating 
particles, where higher concentrations lead to substantially 
larger effects. This happens because the "local" contribution 
is important. 

In summary, the effects of DM decays and/or annihi- 
lations on structure formation are quite small, except in 
some extreme cases (e.g. very high concentration for an- 
nihilations). However, the energy injection from DM de- 
cays/annihilations has important consequences on the frac- 
tion of gas which is retained inside the halo. This fraction 
can be substantially reduced, especially in the smallest halos 
( < 1O 6 M ). 

Finally, we point out that our results are quite different 
from the conclusions of Biermann & Kusenko (2006) and 
Stasielak et al. (2006) 4 , who suggest that sterile neutrino 
decays can favour the formation of first objects. The dis- 
crepancy is likely due to our more complete treatment which 
includes the hydrodynamics of the collapsing structures. In 
fact, our hydro-dynamical treatment allows to describe the 
detailed gas density evolution during the collapse, resulting 
in markedly different temperature and chemical properties 
with respect to those found by a simple one-zone model. 
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